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IMAGING METHOD FOR A MULTI-SLICE SPIRAL CT SCAN WITH 3D 
RECONSTRUCTION, AND A COMPUTED TOMOGRAPHY UNIT FOR 
CARRYING OUT THIS METHOD 

[0001] The present application hereby claims priority under 35 U.S.C. §119 on 
German patent application number DE 102 48 770.7 filed October 18, 2002, the entire 
contents of which are hereby incorporated herein by reference. 

Field of the Invention 

[0002] The invention generally relates to an imaging method for a multi-slice 
spiral CT scan with 3D back projection. An object to be examined may be scanned 
spirally with reference to its absorption behavior by way of a rotating ray bundle 
moving in the direction of the axis of rotation, the measured absorption data may be 
collected by the detector of planar design. The measured absorption data are projected 
onto a virtual detector, preferably at the fulcrum of the rotation, and filtered, and for 
the reconstruction of the absorption value of at least one voxel use is made of those 
measured and filtered data that have been produced by rays which penetrate this at 
least one voxel. The invention also generally relates to a CT unit for carrying out this 
method and having appropriate imaging device. 

[0003] The invention also generally relates to a CT unit for carrying out this 
method and having corresponding imaging device. 

Background of the Invention 

[0004] It is known to make use in computed tomography of multi-slice detectors 
and ray bundles with cone geometry for the purpose of scanning objects to be 
examined, in particular patients, and to carry out a 3D image reconstruction. In this 
process, integral absorption data of rays that penetrate the object to be examined are 
acquired and collected by multirow detectors of planar design, the detectors being 
moved spirally, together with a radiation source, about an object to be examined. 
Subsequently, the measured data are subjected to suitable filtering, and selected 
filtered data are back projected in order to produce tomographic images of the object 
to be examined. In this case, 3D image reconstruction means that for the purpose of 
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back projection individual three-dimensional volume elements (=voxel) are 
considered in each case, and use is made in each case of the data that have been 
produced by rays and which have penetrated this voxel. 

[0005] This mode of procedure is attended by the problem that because the 
collected data, which are firstly projected, for the purpose of simplifying calculation, 
onto a virtual detector - which lies on a plane that intersects the axis of rotation of the 
real detector and radiation source - are present in two dimensions and the filtering is 
intended to take place only in one dimension because of the resultant outlay on 
computation, there is no preferred direction to be seen in which the filtering should 
take place. If an unfavorable filtering direction is selected, it is therefore accompanied 
by losses in image quality. 

SUMMARY OF THE INVENTION 

[0006] It is therefore an object of an embodiment of the invention to find a 3D 
reconstruction method in which the direction of filtering is selected optimally in order 
to achieve high image quality. 

[0007] A consideration of the image of a virtual detector with projections, imaged 
thereon, of segment planes such as are defined below and are described, for example, 
in DE 101 26 638 Al, the entire content of which is hereby incorporated herein by 
reference, shows that the direction of the projections varies over the virtual detector. 
According to an embodiment of the invention, the filtering of the data in the direction 
of these projections is carried out on the virtual detector. If, furthermore, 
consideration is given to the path of the projection of a voxel on the virtual detector 
during rotation, this projected voxel runs through zones of different alignment of the 
segment plane projections. It is particularly favorable in this case to undertake the 
filtering of the data respectively along the alignment of the segment plane projections 
such that different filtering directions are produced as a function of the angle of 
projection. 
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[0008] Considering a 3D reconstruction such as described in DE 101 26 638 Al, 
Ntiit multiply inclined segment planes are fitted into the spiral pitch in a spiral segment 
of length Ls < n+2*3niax (= 180° plus twice the fan angle pmax of the detector), or < 
180° after preceding parallel sorting. The maximum inclination of the segment planes 
Rk (121 <k< Ntiit) is determined in this case from the condition that rays must be 
present for the plane at the two ends of the spiral segment within the measuring field. 
For example, it then holds for Ls = 240° (where pmax = 30°), and for Ls = 180° after 
parallel sorting that N^iit = fs * M/p] , where M is the number of detector rows and p 
the pitch of the spiral. 

[0009] According to an embodiment of the invention, the filtering direction of the 
data is defined by the intersection lines, determined in the spiral segment at the 
spacing Aa = Ls /M, of the segment planes in the virtual detector at the fulcrum. 

[0010] The virtual detector is completely covered if the spiral segment is widened 
on the left and on the right for k = 1 and k = Ntiu, respectively. A separately filtered 
detector image is attained in this way for each segment plane Rk (l<k< Nuu). This 
filter operation is performed in each of the detector images measured by projection, 
the detector images Dk responsible for the back projection for different subsegments 
of the reconstmction interval. 

[0011] A spiral segment Ii of the length Li = [-oCmax, cCmax] where oCmax = M-7i/p, is 

subdivided from the 3D back projection into Ntiu overlapping subsegments 

Ii^ (l<k< Ntiit) of length Ls, their centroids differing from one another by at most Ls 

and being subdivided equidistantly. The segment Ii is fitted centrally into the 

reconstruction segment Ii. This results in the following for the subsegments 

Ir*^ (l<k<Nti,t): 

=I,';l<k<N,, 
=//u{-</ max,-<3f max} 
=// u\amax, a max) 
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such that the projection datum, belonging to an image voxel, in the detector image Dk 
is determined by projection in the reconstruction segment Ir*^, where 
1 <k< Ntiit. 

[0012] It is advantageous to track the image of a voxel in the virtual detector in a 
reconstruction segment Ir in order to illustrate this idea of an embodiment of the 
invention. Considering firstly the spiral segment Ir^ and the projection angle a = - 
ocmax, the detector images Di are responsible for the back projection in this spiral 
segment in terms of projection. All voxels that are imaged in the virtual detector at the 
projection angle a = - oCmax in the selected filtering direction according to an 
embodiment of the invention form a curved surface Hi. 

[0013] The filtering direction at a = - a^ax is given by the image of the segment 
plane Ri centered at the projection angle a = -otmax + Ls/2 and fitted into the spiral 
pitch. The image of the curved surface Hi is not a straight line in the virtual detector 
for oo-ocmax + Ls/2. It has been shown in the publication by H. Turbell, 
P.E. Danielsson, entitled An improved Pl-method for reconstruction from helical 
cone-beam projections, Proc. Int. Conf. On Fully 3D image Reconstmction, Egmond 
and Zee, The Netherlands, June 23-26, 1999, the entire contents of which are hereby 
incorporated herein by reference, that the segment plane Ri constitutes a good 
approximation of the curved surface Hi. Consequently, in all other projections in the 
spiral segment Ir' the image Hi in the virtual detector is grouped around the 
corresponding intersection line of Ri. Thus, the filtering directions of Di that are 
given by the image of R] in the virtual detector are suitable for all the remaining 
projection angles in the subsegment Ir, The above argument is valid in like fashion 
for all the remaining subsegments Ir^ (1 <k< Ntiu). 

[0014] In accordance with this basic idea, the inventors propose to improve the 
imaging method, known per se, for a multi-slice spiral CT scan. It is known that: 

an object to be examined is scanned spirally with reference to its absorption 
behavior i.e. rotating ray bundle moving in the direction of the axis of rotation, and 
the measured absorption data are collected by a detector of planar design. 
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the measured absorption data are projected onto a virtual detector, preferably 
at the fulcrum of the rotation, and filtered, and 

in order to reconstmct the absorption value of at least one voxel, use is made 
of the measured and filtered data that are produced by rays that penetrate this at least 
one voxel. 

[0015] The improvement in the method resides in the particular selection of the 
filtering direction, in that the filtering of the data, used for the reconstruction, in the 
image of the virtual detector is performed in the direction of the projection of spiral 
segments that are imaged thereon and which are produced by the spiral scanning over 
a prescribed angular range. 

[0016] The result of this is that different and in each optimum filtering directions 
occur for the data of different reconstruction segments of a voxel, in which case there 
is the additional advantage that the pitch of the spiral of the selection of the bed feed 
is arbitrary. It is therefore possible to consider overlapping scanning in the z-direction. 

[0017] In a preferred design of the method, the filtering takes place along the 
intersection line of doubly inclined planes in the virtual detector. 

[0018] In accordance with a further preferred design of the method, the prescribed 
angular range is swung - 240'' for a spiral segment of length Ls, the angle 240° 
corresponding to the sum of 180° + 2 x fan angle Pmax that is necessary for the 
reconstruction of an image voxel. 

[0019] Alternatively, it is also proposed to carry out parallel sorting of the rays in 
order to form the virtual detector before filtering, the prescribed angular range for a 
spiral segment of length Ls being 180°. 

[0020] Furthermore, the segment planes formed at least approximately by the 
spiral segments can have a maximum inclination that is limited by the condition that 
rays for the plane must be present at the two ends of the spiral segment within the 
measuring field. 
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[0021] In addition, for the purpose of 3D back projection it is advantageous to 
subdivide a spiral segment Ii of length Li=[-amax, + (W] with Omax = M-7i/p into Nau 
overlapping subsegments Ii^ (1 <k< Nau) of length Ls, whose centroids differ from 
one another by at most Ls, p corresponding to the set pitch, such that the following 
holds for the subsegments Ir*" (1 <k < Ntiu) produced: 
// =//;!< A: <iV,,, 

^/ ^{-(^ max,-6ir max} 

, Ntilt J Ntilt . V \ 

If^ = Ij u\a max, a max/ 
and the projection datum, belonging to an image voxel, in the detector image Dk is 
determined by projection in the reconstruction segment Ir^ (1 <k< Ntiu), a^'max 
representing the maximum angle reached by the ray through the voxel V. The 
projection datum corresponds to the projection value that is found in the detector and 
which corresponds to the image of the image voxel in the detector, the image of the 
point of intersection of the ray from the focus through the image voxel being in the 
detector. 

[0022] If, in addition, consideration is given to the path, of somewhat different 
length, of the rays through the object to be examined as a function of the cone angle 
of the ray bundle used, it emerges as advantageous to weigh the measured absorption 
data as a function of the cone angle of the ray produced in the direction of the axis of 
rotation of the detector and radiation source, preferably with the a cosine of its cone 
angle. 

[0023] The detector of planar design used for the spiral scanning of the object to 
be examined can be provided with a multiplicity of detector elements that are 
arranged matricially in rows and columns. However, there is also the possibility in 
principle of differently fashioned distributions of the detector elements over the 
detector, for example, in a helical arrangement or in an arrangement with detector 
rows of different width that can also, if appropriate, have detectors of different extent 
in the direction of rotation and directional pitch. 
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[0024] In addition to the method according to an embodiment of the invention, the 
inventors also propose a CT unit for scanning an object to be examined. Such a unit 
may be equipped with a ray bundle emanating from at least one focus, and with a 
detector array of planar design and having a multiplicity of distributed detector 
elements for detecting the rays of the ray bundle, the at least one focus moving 
relative to the object to be examined on at least one focus track with the detector array 
opposite. It may further include at least a device for collecting detector data, with 
filtering and 3D back projection being provided. The device for filtering may be 
fashioned in such a way that the method is carried out. The device for filtering can be 
implemented at least partially by programs or program modules. 

[0025] The overall result of the use of the method according to an embodiment of 
the invention and of the inventive CT unit is that an improved image quality is 
attached by the particular selection of the filtering direction. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0026] The present invention will become more fully understood from the detailed 
description of preferred embodiments given hereinbelow and the accompanying 
drawings, which are given by way of illustration only and thus are not limitative of 
the present invention, and wherein: 

Figure I shows multi-slice CT in a perspective illustration of the scanning unit; 

Figure 2 shows a longitudinal section through the multi-slice CT from figure 1 ; 

Figure 3 shows the projection of segment planes in the parallelized 
representation of a physical detector Dp on a virtual detector Dc with a representation 
of three intersection volumes and associated spiral segments; 

Figures 4a - 4d show a representation of the projection of the segment planes on 
virtual detectors with a different projection angle; 
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Figure 5 shows a representation of the profile of projection planes with the aid 
of 16-row detector and a change in the filtering direction; 

Figure 6 shows an illustration of the migration of the projection of a voxel over 
the detector in the case of a progressive projection angle; 

Figure 7 shows a pictorial representation of the juxtaposed reconstruction 
segments of the reprojection of a voxel. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0027] Figures 1 and 2 show a partially perspective illustration of a third 
generation multi-slice CT unit for carrying out the method according to an 
embodiment of the invention. The measuring arrangement (=gantry) denoted by 1 has 
an X-ray source 2 with a beam diaphragm 3 positioned in front of it. This array of 
planar design having a plurality of rows and columns of detector elements forms a 
detector system 5 and is illustrated in section in figure 2 with a beam diaphragm 6 
positioned in front of said detector system and near the detector. For the purpose of 
greater clarity, figure 1 illustrates only eight rows LI to LQ of detector elements 4. 
Detector system 5 can, however, also have any other or preferably greater number of 
rows without departing from the scope of the invention. A different planar 
arrangement of the detectors is also likewise possible. 

[0028] The X-ray source 2 with the beam diaphragm 3 on the one hand, and the 
detector system 5 with the beam diaphragm 6, on the other hand, are fitted on a rotary 
frame 7 opposite one another in such a way that a pyramidal X-ray bundle whose 
edge rays are denoted by eight, which emanates from the X-ray source 2 during 
operation of the CT unit and is stopped down by the adjustable beam diaphragm 3, 
strikes the detector system 5. In this case, the bean diaphragm 6 is set to correspond to 
the cross section of the X-ray bundle, set by means of the beam diaphragm 3, such 
that in the accordance with different operating modes only that region of the detector 
system 5 is exposed which is struck directly by the X-ray bundle. Only eight rows of 
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detector elements 4 are used in figures 1 and 2, while the further rows, indicated by 
dots, are covered by the beam diaphragm 6 and therefore not active. 

[0029] The X-ray bundle has a cone angle cp, which is the aperture angle of the X- 
ray bundle in a plane containing the system axis Z and the focus F. The aperture angle 
of the X-ray bundle is 2/^ax in a plane (an aperture angle) situated at right angles to 
the system axis Z and contained in the focus F. 

[0030] The rotary frame 7 can be set in rotation about a system axis denoted by Z 
by means of a drive device 22. The system axis Z runs parallel to the z-axis of the 
coordinate system which is rectangular in three-dimensions and illustrated in figure 1. 

[0031] The columns of the detector system 5 likewise mn to the direction of the z- 
axis, while the rows run transverse to the system axis Z. 

[0032] In order to be able to bring an object to be examined, for example a patient, 
into the beam path of the X-ray bundle, a support device 9 is provided which can be 
displaced parallel to the system axis Z, that is to say in the direction of the z-axis, 
specifically in such a way as to obtain synchronization between the rotation 
movement of the rotary frame 7 and the translation movement of the support device 
such that the ratio of translation speed to rotation speed is adjustable by selecting a 
desired value of the feed p (= pitch) of the support device 9 per rotation of the rotary 
frame. 

[0033] It is hence possible for a volume of an object to be examined, located on 
the support device 9, to be examined in accordance with a volume scan, it being 
possible to carry out volume scan in a form of a spiral scan such that a multiplicity of 
projections from different projection directions are recorded by way of a measuring 
unit and per revolution of the measuring unit 1, by simultaneous rotation of the 
measuring unit 1 and translation of the support device 9. During the spiral scan, the 
focus F of the X-ray source moves relative to the support device 9 on a spiral track 
denoted by S above the object to be examined. In order to permit a complete 
reconstruction of a CT image per row of detector elements, the spiral scanning must 
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extent in the a direction over at least 7i+2>flinax, but it can also be arbitrary longer 
within the technical limits of the CT unit. 

[0034] However, as a consequence of the circumstance that a plurality of rows of 
detector elements 4 are present, it is also possible for a volume of the object to be 
examined to be examined in the course of a so called tomogram scanning, in the case 
of which no relative movement takes place in the direction of the z-axis between the 
measuring unit 1 and support device 9 (p=0). In the case of the tomogram scanning, 
therefore, the magnitude of the volume examined is determined by the number of the 
active rows of detector elements 4. During a tomogram scan, the focus F is moved on 
a circular focus track or, when the feed is present, a spiral one that lies in a plane 
denoted below as center plane. 

[0035] The tomogram scan can be performed as a partial revolution or in the form 
of a complete revolution, the partial revolution comprising a partial revolution interval 
of at least 7i+2Pniax (a half revolution plus fan aperture angle), thus permitting a 
complete reconstruction of a CT image, while a complete revolution comprises 2n. 

[0036] The measured data which are read out in parallel during the spiral or 
tomogram scanning from the detector elements of each active row of the detector 
system 5 and which correspond to the individual projections P(a,p,q) in fan ray 
geometry are subjected in a data conditioning unit 10 to digital-to-analog conversion, 
serialized and transmitted to an image computer 11. 

[0037] After preprocessing of the measured data in a preprocessing unit 12 of the 
image computer 11, the resulting data stream passes to a tomographic image 
reconstruction unit 13 which uses the measured data to reconstruct tomographic 
images of the desired sections of the object to be examined in a way known per se 
using a 3D back projection method, the data used firsdy being filtered, however. 

[0038] The CT images are put together from pixels (pixel= picture element) 
assembled matricially, by pixels being assigned to the respective image plane, each 
volume being assigned a CT number in Hounsfield units (HU), and it being possible 
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to represent the individual pixels according to the CT number/grey-value scale in a 
grey-scale value corresponding to its respective CT number. In this case, each pixel 
illustrates a voxel (voxel = volume element) of the section of the object to be 
examined illustrated in the CT image. Since measured data are obtained with 
reference to a plurality of sections of the object to be examined owing to the fact that 
the detector system 5 has a plurality of rows and, if appropriate, owing to the spiral 
scanning, 3D data are available, and these are subjected to a 3D back projection 
within the scope of an embodiment of the invention. 

[0039] Available as the end result are 3D image data in the form of a three- 
dimensional matrix, for example with the axes x, y, z, each element of the matrix 
corresponding to a voxel V(x, y. z) and containing the grey-scale value corresponding to 
the associated CT number. Those elements of the three-dimensional matrix that have 
the same x-, y- or z- value then respectively constitute a planar tomographic image of 
the section of the object to be examined that corresponds to the decisive x-, y- or z- 
value. 

[0040] The images reconstructed by the tomographic image reconstruction unit 13 
are displayed on a display unit 16, for example a monitor, connected to the image 
computer 1 1 . 

[0041] The X-ray source 2, for example an X-ray tube, is supplied by a generator 
unit 17 with the required voltages and currents, for example the tube voltage U. In 
order to be able to set these to the respectively required values, the generator unit 17 
is assigned a control unit 18 with keyboard 19 that permits the required settings. 

[0042] The other operation and control of the CT unit is also carried out by way of 
the control unit 18 and the keyboard 19, and this is illustrated by the fact that the 
control unit 18 is connected to the image computer 11. 

[0043] Inter alia, it is possible to set the number of active rows of detector 
elements 4 and therefore the position on the beam diaphragms 3 and 6, for which the 
control unit 18 is connected to adjustment units 20 and 21 assigned to the beam 
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diaphragms 3 and 6. It is also possible to set the rotation time x which the rotary 
frame 7 requires for the full rotation, and this is illustrated by the fact that a drive unit 
22 assigned to the rotary frame 7 is connected to the control unit 18. 

[0044] Although it is also possible in principle to implement the inventive method 
by using fan beam geometry as well, the described CT unit is preferably operated in a 
mode in which the inventive method is implemented using parallel beam geometry. 

[0045] Consequently, the data obtained from spiral or tomogram scanning using 
fan beam geometry, the scanning of the body region of the patient relevant for the 
respective examination are firstly converted in a way known per se into data in terms 
of parallel beam geometry by using a method denoted in general as "rebinning". This 
conversion is based on resorting the data obtained using fan beam geometry in such a 
way that rays are extracted from different projections recorded using fan beam 
geometry and combined to form a projection in parallel beam geometry. In parallel 
beam geometry, data from interval of length n suffice to be able to reconstruct a 
complete image. In order to be able to obtain these data, data must be available in fan 
beam geometry from an interval of length 7i+2Pinax. 

[0046] Figure 3 shows the spiral track S of the measuring arrangement of the spiral 
CT, on which the focus F moves with the physical detector Dp situated opposite, 
conversion having already been carried out using parallel beam geometry in this 
representation. This may also be seen by virtue of the fact that in this image the focus 
F is shown in five juxtaposed positions with three parallel rays in each case. In the 
background, these parallel rays strike the physical detector Dp which, however, is 
already displaying one image of many detector positions of revolving detector 
elements. The detector image of the physical detector Dp is projected, furthermore, 
onto the virtual detector Dv, which adjoins the system axis. Figure 3 also shows three 
exemplary planes of section El to E3, which are imaged on the detector Dv. 

[0047] Considering the virtual detector Dv in the coordinated system solely with 
the representation of the projections PI to P3 that correspond to the planes El to E3 
from figure 3, the result is a representation of the virtual detector Dv such as is shown 
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in figures 4a to 4d. The successive displacement of projections in the figures 4a to 4d 
is produced in this case by considering the planes El to E3 and their projections PI to 
P3 for progressive rotation about the z-axis. 

[0048] The idea of an embodiment of the invention resides also in also carrying 
out in the 3D reconstruction filtering which proceeds in the respective direction of 
projection of the plane considered. As may be seen from figures 4a to 4d, the 
alignment of the projections changes from bottom to top continuously, and the 
filtering direction therefore does so likewise. 

[0049] Consideration of the corresponding representation with a higher number of 
planes and projections results in a representation of the virtual detector as is shown in 
figure 5. The projections P in this case change their direction from bottom to top in a 
typical continuous way. According to the invention, two-dimensional filtering is 
carried out correspondingly for each detector image thus obtained in a fashion from 
left to right, the direction of filtering also being over the height of the detector in 
accordance with the direction in which the projections run. 

[0050] For the purpose of the 3D back projection, which relates to the voxel 
V(x,y,z), the rays that penetrate a voxel in this volume to be examined are considered so 
as to determine the absortivity of a voxel. In the case of continuous revolution above 
the z-axis and projections of this one voxel V(x,y,z) onto the virtual detector Dv, the 
track of the projection of this voxel runs over the multiplicity of the detector images 
considered, as illustrated in figure 6 by the track Bv. Here, the line of section of a 
segment plane in the respective virtual detector is plotted for different directions of 
protection. 32 different directions of projection are considered in this case. The 
inventive method in this case filters the individual detector points in different 
directions f that respectively correspond to the direction of the projections of the 
segment plane which runs through this voxel. It may be pointed out in order to avoid 
misunderstandings that this number is not to be equated with the row number of the 
detector. 
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[0051] The voxel track Bv can be illustrated corresponds to a revolution of the 
measuring arrangement between -Omax'' and H-ctmax^. If this track is subdivided in the 
inventive way into reconstruction subsegments Ir^ to Ir^^'^ these can thus be 
represented as shown in figure 7. It is to be pointed out in this case that the illustrated 
total angle from -amax'' and H-oCmax'' is much larger than the 180° shown and 
corresponds to a plurality of revolutions of the measuring arrangement. With the 
exception of the segments positioned at the edge - segments are arranged 
equidistantly with reference to their centers and centroids, and in a slightly 
overlapping fashion. 

[0052] The back projection of these reconstruction segments onto the voxel V(x,y,z) 
considered then gives the absorption value, assigned to this voxel, inside the object to 
be examined. If this mode of procedure is carried out to all the voxels in the entire 
considered volume to be examined, the result is a volumetric representation of the 
absorption values of the object to be examined. 

[0053] Because if image voxel is imaged in each subsegment of the reconstruction 
region in the detector, each image voxel makes a contribution of the same weight to 
the reconstructured image volume. Consequently, there is no need for the measured 
data to be weighted by voxel. 

[0054] Thus, in sum, an embodiment of the invention uses an imaging method for 
a multi-slice spiral CT with 3D back projection, and a spiral CT in which in order to 
reconstruct the absorption value from at least one voxel, use is made of the measured 
and filtered data that are produced by rays which penetrate this at least one voxel, the 
filtering of the data, used for the reconstruction, in the image of the virtual detector 
being performed in the direction of the projection of spiral segments that are imaged 
thereon and which are produced by spiral scanning over a prescribed angular range. 

[0055] It may be pointed out in addition that the previously described method and 
CT can also be used to represent periodically moving objects, in particular for 
cardiological spiral CT, in which case it is necessary for relevant data to be weighted 
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or selected in a way known per se and correlated with the situation of movement on 
the object to be examined. 

[0056] Exemplary embodiments being thus described, it will be obvious that the 
same may be varied in many ways. Such variations are not to be regarded as a 
departure from the spirit and scope of the present invention, and all such 
modifications as would be obvious to one skilled in the art are intended to be included 
within the scope of the following claims. 
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